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Radiation From a Vertical Antenna Over a Curved 

Stratified Ground 

James R. Wait 

The problem, of a radial electric dipole outside a concentrically stratified spherical 
conductor, such as the earth, is formulated. The solution is facilitated by considering the 
analogous nonuniform transmission line for the radial modes. The general result is then 
transformed to a Watson-type residue or azimuthal mode series, which reduces to the 
well-known result for the homogeneous earth as a special case. Following a method intro- 
duced recently by Bremmer, the residue series is converted to an alternative expansion, 
which is more suitable at short distances. The leading term, of this new expansion corre- 
sponds to the case of the transmitter and receiver over a plane stratified conducting earth. 



1. Introduction 

In an earlier paper by the author [1] l expressions 
were derived for the fields of a vertical electric dipole 
over a plane stratified ground. In a further paper [2], 
the solution was extended to arbitrary antenna 
heights and numerical values of the attenuation 
factor were given. In the latter paper it was men- 
tioned that the effect of the earth's curvature could 
be accounted for. It is the purpose of this paper to 
develop the theory for propagation over a curved 
earth with concentric stratifications. 

The method of solution is a direct extension of 
Watson's [3] method to obtain a solution for the 
electric dipole radiating over a homogeneous sphere. 
It therefore does not seem necessary to outline the 
steps in the derivation in detail. A particular novel 
feature of the formulation, however, is that the 
boundary conditions at the surface of the sphere are 
specified by a surface impedance. The final expres- 
sions for the fields are then also applicable to propa- 
gation over other modified surfaces, such as a sphere 
with azimuthal periodic corrugations. Using a 
method suggested recently by H. Bremmer, 2 an 
expression for the field is then obtained that is very 
suitable for computation when the receiver is near 
the optical horizon. Finally, some numerical results 
are presented for various frequencies employing 
typical values of the ground constants over both a 
homogeneous and stratified earth. 

2. Formal Solution 

The source of the field is considered to be an elec- 
tric current element Ids oriented in the radial direc- 
tion to the spherical earth of radius a x . Choosing a 
spherical coordinate system (r, #,</>), the surface of the 
earth is then defined by r=a Xl and the dipole or 
current element is located at r=b and 0=0. The 
fields can now be expressed in terms of a Hertz 



1 Figures in brackets indicate the literature references at the end of this paper 

2 Personal communication. 



vector, which has only a radial component U, as 
follows: 

where &=2x/free-space wavelength, and e is the 
dielectric constant of free space, ( = 8.854 X10~ 12 ), 
in mks units. A time factor exp(ico^) has been 
implied . The function U satisfies the equation 



2wr z sin 6 



(2) 



where the 5's are the Dirac delta or impulse function. 
The factor 2irr 2 sin 6 is the Jacobian of the trans- 
formation from rectangular to spherical coordinates. 
The constant is chosen so that U has the proper 
singularity at the dipole, that is 



U- 



AiricjoeR 



Ids for i?->0, 



(3) 



where R=[r 2 -\-b 2 -2br cos 6} 1/2 , 6=8.854 X 10" 12 , and 
tnerefore C= (i/ we) Ids. 

The field U is now written as the sum of two parts 
U e -{-U s , where U e has the proper dipole singularity 
at R=0, and U s is finite at that point. As U s 
is a solution of the homogeneous wave equation, it 
can be written in the form 

U s ^-p^(2n+l)AM ) (kr)P n (cos d), (4) 

where h™ (kr) is the spherical Hankel function of the 
second kind, which assures outgoing waves at 
infinity, and Pn (cos 6) is the Legendre function. 
The corresponding representation for U e is [4} 
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U e = l i±j:(2n+l)h™(kr)h™(kb)P n (cosd) for r<b 

(5) 
= % ^j:(2n+l)h™(kr)hy(kb)P n (cosd) for r>b 

The coefficients A n are now found from the boundary 
condition, that 



Ee=—ZH (f) at r=a u 
which can be rewritten 



L^ or _r= a i 



(6) 



(7) 



In other words, it is assumed that the surface of the 
earth exhibits the property of surface impedance. 
Z is taken to be equal to the ratio of the tangential 
electric and magnetic fields for a vertically polarized 
plane wave at grazing incidence on a plane stratified 
earth. This step in the analysis leads to a great 
simplification, and it is justified in the appendix. 
It then readily follows that 



A n * 



Ai u (tei) 



hf{k(H) 






h?{kb), (8) 



where A=ecoZ/& = Z/1207r and where h^\x) is the 
spherical Hankel function of the first kind. The 
total field is then of the form 



U=J2(2n+l)f(n)P n (cose). 

71 = 



(9) 



Following the process developed by Watson [3], the 
summation is transformed into the following contour 
integral: 

where the contour Ci-\-c 2 encloses the positive real 
axis, as illustrated in figure 1. Noting that the 
poles of the integrand are located at n=l/2, 3/2, 5/2, 
. . ., etc., it can be readily verified by the theorem of 
residues that this integral is equivalent to eq (9). 
Now, since /(ti— 1/2) is an even function of n, the part 
of the contour Ci above the real axis can be replaced 
by the contour c[, which is located just below the 
negative real axis. The contour cJ+^2 is now entirely 
equivalent to i, a straight line running along just 
below the real axis. Replacing n— 1/2 by v, the 
contour representation for U takes the form 



Imag. n 



p "'Ji 



L sin vir 



j(v)P p [cos(T-6)]dp. (11) 



It is to be noted that this manipulation of the con- 
tours is simplified because /(n— 1/2) is an even 



r 



.__±j + .• 



v 
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Figure 1. Complex n plane showing the contours and the zeros 
of cos nw and zeros of the denominator of equation (8) . 

function of n as a consequence of the method of 
formulation. In the usual treatment for the homo- 
geneous sphere, j(n— 1/2) is not an even function 
of n, and the deformation of the contour becomes 
more intricate. 

The next step in the analysis is to close L by an 
infinite semicircle in the negative half -plane. The 
contribution from this portion of the contour vanishes 
as the radius of the semicircle approaches infinity. 
The reasoning for this fact follows directly from 
Watson's argument for the homogeneous sphere. 
The value of the integral for U along the contour L 
is now equal to the sum of the residues of the inte- 
grand evaluated at the poles v s oif(v) located in the 
lower half -plane. It then follows that U is propor- 
tional to 



s = 



(v,+\/2)h%(kb)h%(kr)P, t [coa(T-0)] 



~d{v)~\ 



(12) 



where the function M(v) is defined by 

M(v)={ 1 -~logxh^(x)~if\ (13) 

and the poles v s are the solutions of M(v) = 0. Again, 
as a result of the formulation, the equation here for 
the determination of the roots is relatively simple. 
Making the usual approximation [1] that h™ (x) can 
be represented by Hankel functions of order 1/3, eq 
(13) can be replaced by 

/3 ^[(i/ 3 )(-2x 8 n = _ 2 



where 



8= — i 



1207 
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(AfflO'/'Z 



and r s = 



v s — ka-i 



(ka,y 



(15) 



Tor a homogeneous earth whose propagation con- 
stant is 71, it follows thai 



7j L 7? J 



with jo=ik, and therefore 



Tl/75 



(te,) 1/3 [M/T2)-i] 



1/2 7 



(16) 



(17) 



which is identical to the value given by Bremmer [4]. 
For a two-layer ground whose upper stratum of 
thickness I has propagation constant y u the surface 
impedance Z is given by (see appendix) 



where 
and 



Z=120^ri-^TV, 



<2=tanh[P+(7?-7§) 1/2 Z], 



■ tonh "'[| W1 



(18) 
(19) 
(20) 



with 72 as the propagation constant of the lower 
medium. Tin 4 correction factor Q, which approaches 
unity for I— >co } can be evaluated with the aid of 
charts of hyperbolic functions of complex argument. 
The expression for Q for a stratified ground for any 
number of layers has been given previously along 
with some numerical values for special cases [5]. 

Equation (12) for U can now be considerably 
simplified by replacing the Legendre function by 
its leading term in its asymptotic expansion and the 
functions h™(Icb) and A?(ir) by their Hankel 
approximations to lead to 



U=2UoT/-2iriX J2 MhOMk) 7T- 

* = * &Ti 



nth 
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(21) 



(22) 



and where 

TT Idse- ika i° _ 
c/o=i— = *> hi=r—a u 



47r/coeai(9 



h 2 =b~a u X=(kaO u % 



and 



X^ikatf'^hija^ 2 for 7 = 1,2. 



The preceding equations are identical in form to 
that obtained for the homogeneous earth as given 
by Watson [3], Bremmer [4], Friedman [6], and 
others. It is important to note that the quantities 
p, 8, and the roots t s are dependent on the electric 
constants and the nature of the stratification of the 
earth. Bremmer [4] has given very adequate for- 



mulas for r s in terms of powers of 8. These can be 
used directly for the stratified earth. They can also 
be used for a sphere with a corrugated surface if the 
appropriate value of the surface impedance is 
employed [7]. 

3. Modified Flat Earth Formula 

The so-called residue series for U could be used for 
calculations for propagation over a curved earth. 
The series, however, converges poorly for distances 
near the optical horizon. It would be desirable to 
transform the residue series formula to a new type 
of expansion where the first term corresponds to 
the radiation of a dipole over a plane stratified 
earth. Succeeding terms would then be preferably 
in proportion to inverse powers of ka x . In the 
limiting case when ka x tends to infinity, the expres- 
sion U should correspond to the situation treated 
previously. A method of obtaining expansion for- 
mula of this type is mentioned briefly by Bremmer in 
his book [4]. Very recently he has described to the 
author an alternative procedure, which he illustrated 
for the case when the transmitter and receiver are 
located on the surface of a homogeneous spherical 
earth. His method will be employed here in the 
case when the transmitter and receiver are not both 
on the earth. Furthermore, as will be shown, the 
method also is applicable to stratified and corrugated 
surfaces. The first step is to express the field in 
terms of a contour integral as follows: 



JUJA 

2p l/2 U, 



2*ij e - 



c + i<° _l/2 



[ l+ e ^'Y /2 A(p) 



(23) 



where p—iX/2S 2 and where" \A(p) and f(p,h) are 
chosen as follows: 



A(P)- 






(24) 



and 



■*?+(p/«TF fffflK-33+p/OT 
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(25) 



for i=l, 2; and c is some positive constant. 

It can be noted that the poles p s of the integrand 
are determined by the solution of 



H?i 



iZ 



m 



+e , *"pl' 2 fl$ 



[f]=< 



(26) 



li p s is replaced by 28 2 r s e +iir this equation is identical 
to eq (14) for the determination of the roots t 8 . It 
can be readily verified that the sum of the residues 
evaluated at the poles p s leads back to eq (21). 

Recognizing that the right-hand side of eq (23) is 
in the form of an inverse Laplace Transform, it can 



239 



be written in Heaviside's operational notation [8] 
as follows: 



iU(p)\ 



=*Xp), 



where 



and 



F(j>) 



2p l/2 U ~ 
- i + e^Y^Aip) 



(27) 



(28) 



flforp>0 
[oforp<0. 



The technique is now to expand F{p) in an asymp- 
totic expansion in powers of 1/p and then to invert 
each term to obtain series expansion for p in powers 
of <5 3 . It is convenient for the moment to consider 
A 2 =0 so that/(p,A 2 ) = l. Each of the Hankel func- 
tions is now expanded in an asymptotic expansion, 
with due regard being paid to the phase of arguments. 
Note that 

3 

+ tt>- arg. p s — 3 arg. 8— 3tt> — 2t, 

which suggests replacing the arguments z by ze~ ziir 
and then employing the relation 



m 2 \ze- u «)^~ %1 ^ 2 ™ g»(g). 



sm ttv 



(29) 



On making this conversion, the validity of the 
asymptotic expansions is established. It then follows 
that 



F(p)= 



TT^pfipA) 
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/(p j A 1 )=e- (2 '-^"' 
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h} (31 > 
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plus terms in ^ etc., with 



x? lTV 



and 






To further simplify the above expression for j(p,h) T 
the factor (l + .Y 2 6 2 /^) 3/2 is expanded in a binominal 
series, and terms containing Xj, X\, etc. are neg- 
lected. After some algebraic manipulation, it follows 
that 

f( V ,h0^e «»' |^i-_ — —J (32) 

plus terms containing a 2 5 6 , a 3 <5 9 , etc., where a=khiA. 
When this formula for f{p y hi) is substituted into eq 
(30), it is readily shown that 



F(p)__ iper aV 



IT 



i-ip 1/2 ^ I2(i~ip l/2 y 2 2(i-i v v 2 )y 

5«p? I ! + a 4 

Ls p z/2 {\~ip 112 ) 2 4p(i-ip l/2 Y^4p(i-ip l/2 y n sp 3 



-N 



plus terms containing 5 9 , <5 12 , etc. 

The final step in the analysis is to find the cor- 
responding functions of p for each term on the right- 
hand side of the above equation. The necessary 
operational pairs can be derived from the basic re- 
lation [9], 

_ e -apV2 

Af(g.P)l= 1 j_^w -l (34) 

where 



'H\-iv ll2 )\ {,i6) 



■ + 0P 



with 



M( 9 , P ) = erfc (^) -e^) erfc (^+f ) 



"f/x. 



«*> 1/2 r d 
F 7 T = L^' lf( ' / ' p)+A% ' 



P) 1 (35) 



For example, 

(1 
so that 

-iT 1/2 (l + 2 P -ia)e- l "- t " ) erfc [^ u - 2 +ip V2 ~\ \l (36) 

The other pairs can be derived in a similar manner 
by further differentiations with respect to the param- 
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eter g and then setting g 

collecting; like terms is 



The final result after 



^ =e -a'/4p_; (rp) l/2 e 
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(37) 



plus terras containing <5 6 , <5 9 , etc. 

The coefficient of the <5 6 term is not written out here 
as it is quite cumbersome and in itself would not be 
suitable for computation. It is constructive at this 
stage to express U(p) in terms of the parameter w 
defined by 

w =o( 1+ ^) 2 ' (38) 

whence 



U(p} 
2U 



^j£-ikh\/2a x d J 



l~i(T P ) 1/2 e- w erjc(iw 1/2 ) 

e~ w erfc(iw l/ 



+5 3 r p -i(7r P ) 1 / 2 ^^e-^r/c(w 1/2 )l\ 



+^i(^ <1+i « l4) „ / ,( i «i)"«] (m 

For the homogeneous earth, the factor A= 
(to/ti)(7i/7o— l) 1/2 7 an d if a>i the radius tends to 
infinity or <5=0, the remaining first term corresponds 
to the well-known formula for the attenuation factor 
of a dipole over a flat earth [10]. p is then the 
numerical distance of Sommerfeld [11]. It is inter- 
esting to note that the form of the first term is 
identical to the result obtained by Hufford [12], 
using an integral-equation formulation. When the 
ground is stratified such that A=euZ/k, where Z 
is the surface impedance, the first term corresponds 
to the result obtained previously for the dipole over 
the plane stratified earth [1, 2]. 

It now appears that, for small heights such that 
kh\j2aB is small compared to one, the height gain 
function is simply l-\-ia, which is a common factor 
of at least the first three terms in the expansion in 
powers of <5 3 . To this approximation the final result 
is expressed conveniently as 

^ = IF(p) = ^^o(p)^[1-^p) 1/2 -(1 + 2p)^o(p)] 

+ 5 6 [l-i(7rp)^ 

(40) 
where 

F (p)=l-i(irpy /2 erfc(ip l/2 ) 

#~(1 + ia)(l + ip) with p=kh 2 A. 



At distances from the antenna greater than a few 
wavelengths the quantity U(p) is proportional to 
the vertical component of the electric field, and cor- 
respondingly 2f7 is proportional to the field of an 
identical source over a flal perfectly conducting 
earth. The field strength E, in millivolts per meter, 
at a distance D km , in kilometers, is (hen given by 



E- 



300 

Is km 



\W( P )\ 



(41) 



for a dipole transmitter whose strength is such that 
it would radiate 1 kw over a perfectly conducting 
flat ground. 

To illustrate the nature of the results, calculations 
are carried out for E as a function of D km over a 
homogeneous spherical earth of radius 6,373 km 
at frequencies of 30, 150, and 750 Mc. These curves 
are illustrated in figure 2, where it is indicated that 
the relative dielectric constant of the ground is 4 
and the conductivity is 10~ 2 mho/m. The flat- 
earth formula (i. e., <5->0) is designated as the 
zero order or Oth approximation. The results for 
including the first and second curvature corrections 
are designated as first- and second-order approxima- 
tion, respectively. Finally, the field strengths ob- 
tained from the residue series are shown for compari- 
son and designated by "Res." In this latter case 
many terms of the series were required, and it is 
reassuring to notice that curvature corrected flat- 
earth formulas merge with the residue-series formula. 
The results of these calculations would lend confi- 
dence to the use of eq (40) for stratified media if the 
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Figure 2. Field-strength versus distance curves for a trans- 
mitter that radiates 1 kilowatt. 

The flat-earth and curvature corrected curves are shown along with the residue 
series-calculated curve. 



In eqn. (40), the definition of F ( p) is 

1 °i 

F Q (p) = 1 - i (tt P ) 2 e" p erfc (ip 2 ) . 
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Figure 3. Field-strength versus distance curves for various 
receiving-antenna heights. 



appropriate value of the numerical distance p and 
curvature parameter 8 are employed. As a further 
check on the curvature corrected formula used in 
conjunction with the height-gain function, calcula- 
tions are carried out for various receiver heights 
A 2 with the transmitter height hi equal to zero. 
These results are shown in figure 3 for a frequency 
of 43 Mc, relative dielectric constant of 4, con- 
ductivity of 10~ 2 , and h 2 =0, 5, 10, 20, 50, and 100 m. 
The solid curves correspond to the second-order 
curvature corrected flat-earth formula, whereas the 
dashed curves are based on the cumbersome but more 
exact residue series [4]. As a matter of interest, 
the line corresponding to the optical horizon is shown. 
It is apparent that the agreement between the two 
methods of calculation is very good at moderate 
ranges. No doubt a third-order correction involving 
a term proportional to <5 9 would improve the dis- 
crepancy at greater ranges. It is doubtful, how- 
ever, if calculations of higher order correction terms 
are justified because the residue series becomes more 
convenient at larger ranges. 

The computations for stratified ground are some- 
what more involved, and, furthermore, even for a 
two-layer ground, additional parameters to be con- 
sidered are conductivity a 2 and dielectric constant e 2 
of lower stratum, and thickness I of upper stratum. 
As an example, the lower medium is taken to be 25 
times better conducting than the upper stratum. 
For convenience, the dielectric constants are also 
taken to be in the same ratio. The second-order 
curvature corrected results for a frequency of 5 Mc 
is shown in figure 4, where the case l=oo corresponds 
to a homogeneous ground of electric constants o- x and 
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Figure 4. Example of the field-strength-distance curves for a 
two-layer stratified ground for several values of the thickness 
of the upper stratum. 
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Figure 5. Field strength at a distance of J+O kilometers on a 
two-layer stratified ground as a function of the thickness of 
the upper stratum. 

ei. The curve marked Res. corresponds to the residue 
series calculation for the homogeneous earth. If the 
thickness of the upper stratum is greater than about 
2.5 m, the field corresponds quite closely to that 
expected over a homogeneous ground. To show in 
a little more detail, the nature of the dependence on 
stratification, the field at 40 km is shown plotted in 
figure 5 as a function of I in meters. Both the 
flat-earth and the second-order curvature corrected 
results are included. 
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4. Concluding Remarks 

It is apparent that the theory of Watson [3], 
Bremmer [4], and others, for propagation of radio 
waves over a spherical homogeneous earth can be 
generalized to propagation over a concentrically 
stratified sphere. The derivation in this paper is 
very straightforward as a result of choosing an 
approximate boundary condition, which specifies 
the surface impedance at the surface of the sphere. 
Although the theory was developed for a homo- 
geneous atmosphere, there is little reason to doubt 
that a similar analysis can be carried out for an 
inhomogeneous atmosphere following the develop- 
ment of Friedman [6]. For the case of normal 
refraction (absence of ducting), the results in the 
present paper can be employed directly by admitting 
the concept of the effective earth radius [4], where a { 
is replaced by a modified value a\. 

While the attention here has been mainly confined 
to propagation of radio waves over a stratified sphere, 
the results are applicable to the study of surface 
waves on corrugated spherical surfaces of the type 
that can be characterized by a surface impedance. 

5. Appendix 1. Surface Impedance of a 
Stratified Conducting Sphere 

The surface impedance at the surface of a sphere 
of radius a x consisting of a concentric core of radius 
a 2 will be derived. The propagation constants of 
the outer shell and the core are Ti and 7 2 , respec- 
tively, and the intrinsic impedances are rji and rj 2 , 
respectively. Under the assumption that the source 
is a vertical antenna the field can be represented as 
a superposition of transverse magnetic (TM) spheri- 
cal waves of order n. The surface impedance Z n for 
such a TM wave is defined by 



Z'=[-Et/Hl] r . 



(42) 



The medium between the limits' r=a x and a 2 can now 
be regarded as a nonuniform transmission line [13] 
of length l=a x —a 2 . The characteristic impedance 
of the line looking inward is 






(43) 



where /„ is Schelkunoff's notation for modified 
spherical Bessel functions of order n, and the prime 
indicates a differentiation with respect to the argu- 
ment y { r. Similarly, the characteristic impedance 
looking outward is 



n¥^=N^ 



(44) 



i£n(T 2 r) 
where K n is the modified spherical Hankel function 



of order n. The line is considered to be terminated 
by an impedance 






(45) 



From the theory of nonuniform transmission lines it 
follows that 

where g e and g h are reflection coefficients at r=a 2 
for the El and Hi field components and are given by 



_ 1/P( g2 )-1/M( a2 ) 
ie 1/P(a 2 )+W(a 2 ) 



and 
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P(a 2 )-M(a 2 ) 
= P(a 2 )+N(a 2 ) 



(47) 



(48) 



The quantities A and B are transmission factors 
given by 
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(49) 



Under the restriction that the thickness of the shell 
(l=a 2 —ai) is small compared to a u the above for- 
mulas can be greatly simplified. For example, 
noting that I n (x) satisfies the equation 



d 2 I n (x) 



dx 2 



[l+^]/>), (50) 



= H 



it readily follows thai the function M(x) satisfies 



[i+^]4*=^- 



dx 



A good approximate solution is 

n(n+l)l 1/2 



M( ^{ 1+ wl 



(51) 



(52) 



because dM/dx is small when |7ir|»l. Now the 
values of n are not known precisely, but the impor- 
tant ones correspond to the roots v s , which are close 
to the value ka x . It therefore follows that 

M(^,{l + *M]%[l-|J'', (53) 
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since k = iy , ka{^>\ and di/rc^l in the range 
a 2 <r<a 1 . With similar reasoning it follows that 



r y 2 "1 1/2 
NiJir^nA 1 — I for a<r<a 



and 



P( m )~, 2 [i_|] 



2-11/2 



The transmission factors are then written 
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(56) 



with l=ai— a 2 . The surface impedance then takes 
the form 

arvi— 
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tanh [(Tf-T 2 ,) 172 /] 



(57) 



It is to be noted that Z n is precisely the ratio Z of 
the tangential electric and magnetic fields for a 
vertically polarized wave at grazing incidence on a 
two-layer stratified ground [5]. It, therefore, seems 
justified to employ the approximate boundary eq (7) 
at least for application to propagation of radio waves 
over the surface of a stratified ground with the normal 
earth curvature. 

6. Appendix 2. Alternative Approach to 
the Height-Gain Function 

It is possible to study the effect of raising the re- 
ceiver and/or the transmitter by starting with a 
more accurate representation for the height-gain 
function: 






where 



v=ka-\-(ka] 



1/3 



L 2S'J 



(58) 



with r=h t -\-a; f(p,h t ) is now expanded in a 
MacLaurin series as follows: 



J(p,h)=l+khi 



dx 
.~h 



v \X) Jx = tea 

2 L h?\x) 



+ 



(59) 



Now employing the relation 
and the differential equation 
it follows that 



dx 2 x dx 



j(p f ht)^+(i^^+^\ 



h t \ iAk 2 h\ , k 2 h 2 A 



ka 



(60) 



(61) 



(62) 



It then follows that a more representative height- 
gain function G' ' (p) in place of the function G is 
given by the operator 



G'(p)~l + ikh t A-i 



k 2 h\A , k 2 h\A d 



ka 



2 di 



+ 



(63) 



The field E(p,h) is then related to the field E(p,0) 
on the ground by E{p,h) = G[(p)E(pfi). The third 
term of the expansion is negligible when kh<^(ka) 1/d 
because |A| is somewhat less than unity and the 
fourth term is negligible when 

miD<i. 

The function G'(p) to this approximation is then 
identical to the height-gain function G obtained 
earlier. 
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